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^^ ■ In this paper, we examine the mean field electronic structm'e of the /^ 

j^ . /^ Anderson lattice model in a slave boson approximation, which should be 

useful in understanding the physics of correlated metals with more than one f 
^^ , electron per site such as uranium-based heavy fermion superconductors. We 

P^ I find that the multiplet structure of the /^ ion acts to quench the crystal 

O ■ field splitting in the quasiparticle electronic structure. This is consistent with 

experimental observations in such metals as UPt^. 
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Localized electron systems exhibit multiplets of definite L, S, and J due to the symmetry 
of the Coulomb interaction about an atomic site. The degeneracies of these multiplets 
are not properly described by standard electronic structure techniques used in solids which 
assume the filling of single particle levels using Fermi-Dirac statistics. Despite this, a Landau 
quasiparticle approach appears to be successful even in metals which exhibit high energy 
multiplet excitations. A good example of this is the heavy fermion superconductor UPt^. 
This metal exhibits canonical Fermi liquid behavior at low temperatures with a Fermi surface 
in reasonable agreement with band calculations.!!! On the other hand, high energy neutron 
scattering data reveal an excitation at about 0.4 eV,a close to the 0.5 eV transition seen in the 
free ion, known to be a multiplet excitation between the '^H^ ground state and the ^F2 excited 
state. Therefore, UPt^ is quite unusual in that its low energy excitations exhibit single 
particle Kramers degeneracy whereas its high energy ones exhibit multi-particle Coulomb 
degeneracy. Our purpose here is to understand what infiuence these p multiplet correlations 
have on the low energy quasiparticle excitations. We do this by studying the f^—p Anderson 
lattice model at the mean field level within a slave boson approximation. Since this study 
deals with the relation of single particle statistics to multi-particle statistics, it should be of 
relevance to other problems as well. 

The /^ — /^ Anderson model has been investigated by a number of authors. Early work 
studied the single impurity model either from a variational approacfm or by use of the 1/N 
expansion!^ where N is the degeneracy of the ion. The lattice case was in turn looked at 
by various authors within a slave boson approximation,!^}!^ but the effects of multiplets were 
ignored. Multiplet notation was introduced in what we will denote the mn basis by Newns 
and Read,!^cl but no calculations were carried out explicitly including multiplet effects (the 
mn state has one f electron in level m and the other in level n). This work was extended by 
Dorin and Schlottmann§ who dealt with the general P - /-^ problem in the mn basis and 
also included Gutzwiller corrections to the mean field approximation, which are necessary 
when looking for possible magnetic instabilities.!^^ Our treatment in the mn basis will closely 
parallel this work. In turn, though, we will develop the formalism in the JJz basis, where 
J is the total angular momentum of the /^ state. This basis is diagonal in the Coulomb 
interaction and contains the proper multiplet degeneracy of the /^ ion. 

The ionic part of the Hamiltonian in the mn basis is 

y^ Em\m >< m\ + ^^ Emnl^TT^n >< 'fnn\ (1) 

jm jmn 

where Em is the energy of one f electron in level m, Emn is the energy of two f electrons 
in levels m and n, and j is the site index. Note that the state \mn > is antisymmetric in 
the two indices. We will now assume that we are in the strong spin-orbit coupling limit (as 
appropriate for uranium) so that the /^ state has a J of 5/2 (6 states) and the /^ state has 
a J of 0, 2, or 4 (15 states). 

The conduction electron part of the Hamiltonian is 

/ ^ km %iYi ^"^ ^ ' 

km 

where m is assumed to run over the same 6 indices as the /^ state. 
The hybridization part of the Hamiltonian is 



j:ykm^"''''\^n><n\c^,^ + H.C. (3) 

jkmn 

where H.C. denotes the Hermitian conjugate. In this paper we replace Vj:^ by a constant, 
V. 

To solve this Hamiltonian, we employ the slave boson approximation. We associate a 
boson pjm at each site j for each /^ state in level m and a boson djmn at each site j for 
each p state in level mn. For the mean field approximation we consider, these bosons are 
replaced by c numbers assumed to be site independent. The ionic part of the Hamiltonian 
per site can now be written as 

E^™Pm + E^-"^L (4) 

m mn 

We must also satisfy a completeness relation 

EPm + E^mn = l (5) 

m mn 

We associate a Lagrange multiplier A^ with this constraint. The other constraints involve 
the f count in each channel m at each site 

i^m = JmJm = Pm ' 2-^ "'mn V"J 



n 



with which we associate Lagrange multipliers A^- 
Finally, the hybridization term now becomes 



E yfy^m^^ + H.C. (7) 

where 



km 
km 



Z< 



m 



E '^rnnVn (8) 



n 



and /^^ is the Fourier transform of jjra- The Gutzwiller generalization would involve re- 
placing Zra by£3 

(9) 



'-m I 

We do not use such a renormalization here, but it can be easily included. 

We now have to determine the following parameters: 6 pm, 15 dmn, A\ and 6 A^. These 
are determined by setting the derivative of the Hamiltonian with respect to each of them 
equal to zero. The final parameter is the chemical potential, fi, determined by the total 
number of electrons. 

The quasiparticle bands which come out of this model are the standard ones 

Et = ^'^^ ± ^\/(^'" - '>^J' + 4-^^' (10) 



The following mean field relations are necessary to solve these equations (we assume a 
constant conduction electron density of states, NqjB 

£ < 4„<=*™ >= '^»('' - c) (11) 






E ^e < 4c,„ >= NoZ^V ln(^^) (13) 

g t„ Am 

which we denote by /„, Jm, and i^^, respectively where e^ is the bottom of the E^ band. 
These equations as written assume only the E~ bands are occupied for notational simphcity, 
but in actual calculations, we used the full expressions which can include contributions from 
the £"+ bands depending on the value of the chemical potential, /i. We note that the f count, 
Hf, is simply the sum over m of the Jm- 

The mn basis is not diagonal in the Coulomb interaction, so to treat the effect of multi- 
plets properly, we convert the above formalism to the JJz basis. We note that 

\JJz>=Y.UTj:\mn> (14) 

mn 

where the U are Clebsch-Gordon coefficients.t3 These act to define a new boson field for the 
P state which we denote as djj^. 

The ionic part of the Hamiltonian per site is now 

Hionic = Yl ^raVm + Yl ^J-h^JJ, (l^) 



with the completeness relation 



Ep^^ + E4.. = 1 (16) 



To convert Eqs. 6 and 8, though, one must remember that d is actually an operator and 
that dmn and djj^ represents its expectation value in each basis. Thus, taking advantage of 
the assumed diagonality in the JJz basis, we have 

dmn = j:UTZdjJ^UTf^ (17) 

JJz 



dLn = Y.u7FAjUri (18) 

JJz 

These two equations are then inserted into Eqs. 6 and 8 to obtain the values of hm and Zm 
in the J J, basis. 



Solving the full problem above is quite difficult given the large number of variational 
parameters. We will therefore consider a simple model which illustrates the relevant points. 
In this model, we only keep one /^ configuration, the doublet J, Jz = 4, ±4. Since the /^ 
state I ± 1/2 > does not couple to this, we elect to keep only the | ± 3/2 > and | ± 5/2 > 
states in the /^ and conduction electron manifolds, thus further reducing the number of 
variational parameters. Since we are considering an unpolarized solution, one has only the 
following equations to solve 

2(p^/2 + P5/2 + 0-1 = (19) 

pI/2 + dl, - J5/2 = (20) 

pI/2 + 4 - h/2 = (21) 

(-E44 + A — A5/2 — A3/2)c?44 
+ V^(^5/2P3/2 + i^3/2P5/2) = (22) 

(^5/2 + A^ - A5/2)p5/2 + VK^/2d^^ = (23) 

(E3/2 + A^ - A3/2)P3/2 + VK^/2dii = (24) 

2( J3/2 + J5/2 + /3/2 + h/2) -Nt = (25) 

with 2:5/2 = (^44^3/2, 2:3/2 = (^44^5/2, and Nt the total number of electrons. We solve these 
seven equations by squaring them all and adding them together. The resulting master 
equation is a function of seven variables and is minimized using Powell's method.ta 

For the results presented, we took a conduction band width of 2 eV, a hybridization 
potential, V, of 0.7 eV, an E5/2 of -0.5 eV, an E3/2 of -0.4 eV, and a total electron count 
of 3.5. £'44 was varied in steps of 0.05 eV. Stable solutions were found for ,^44 0.13 eV and 
lower. For higher values of ii^44, the minimiziation failed to converge from a typical input. 
A solution at 0.15 eV was generated by trying a large number of input parameters. For 0.20 
eV, the lowest minimum was found for an n/ of 1. For low values of E44, one asymptotically 
approaches an nj of 2, so the runs were terminated at -1.50 eV. 

The results are presented in the three figures. Fig. 1 shows a plot of the renormalized 
crystal field splitting, AA = A3/2 — A5/2, versus Uf. The data for nj > 1.2 appears to follow 
a straight line which goes to zero at n/ = 2 and interpolates to half the bare crystal field 
splitting at Uf = 1. Below an Uf of 1.2, the chemical potential jumps from the lower to the 
upper 5/2 band. After this, the crystal field splitting begins to approach the unrenormalized 
value of 0.1 eV (the lowest minimum for E44 of 0.20 eV had an nj of 1 and an unrenormalized 
crystal field splitting). We contrast this behavior with the /° — /^ Anderson lattice model, 
where no renormalization of the crystal field splitting occurs. In Fig. 2, a plot of A5/2 — /i is 
shown, which should be an estimate of the Kondo temperature for this model (the absolute 
value is shown, since this quantity changes sign near an nj of 1.2). This vanishes linearly 
with Uf as Hf approaches 2 but vanishes faster than this when Uf approaches 1, with the 
latter behavior beginning when Uf drops below 1.2. This same type of asymmetry was seen 
in earlier impurity calculations.Q In Fig. 3, we show the f occupation in each channel as 
a function of nj. Again, a change in behavior occurs when nf drops below 1.2 with very 
different limiting behavior as nj approaches 1 or 2. We suggest, then, that the jump of the 
chemical potential is necessary for the solution to cross over from the behavior characteristic 
of Uf near 2 to obtain the unrenormalized behavior as Uf approaches 1. The value of Uf at 
which this occurs is model dependent. 



We now present a derivation of the linear behavior of Fig. 1. Subtracting Eq. 24 from 
Eq. 23 and utihzing the defintion of K^n in Eq. 13, we have 

AA = AAo + N.VMl, h \'^'~^%'^"' ] (26) 

^5/2 - ^^ea - A3/2 

where AA = A5/2 — A3/2 and AAq = -E5/2 — -^3/2- As nj approaches 2, we can set the second 
fraction in the argument of the logarithm to 1, and expanding the logarithm, we obtain 

A A ^ AAo - N^V^dl^-^^ (27) 

where \av = |(^5/2 + ^3/2)- This approximation has been checked numerically and works 
well for Hf near 2. To get an expression for A^^, — /i, we analyze the definition of Jm in Eq. 
12. Since Zm is vanishing as Uf approaches 2, this forces A^ to approach /i, and so we can 
ignore the first term in parenthesis in Eq. 12. Using Eqs. 20 and 21, we obtain 

2 
A..-/x^iVo\/^4 /7,, (28) 

where pl^ = 1{pI/2+pI/2)- From Eqs. 19-21, 4pl^ = 2 - n/, so 

A,, - ^ ~ NoVdj,'^—^ (29) 

Substituting this into Eq. 27, we obtain 

AA^^(2-n,) (30) 

which is the desired relationship. 

Fig. 1 is the main result of this paper. We note that except for the special cases of 
the J = 4, Jz = ±4, ±3 states, all /^ states involve combinations of states in all three of 
the /^ subbands. Because of this, in general, all crystal field splittings in the quasiparticle 
electronic structure should be expected to vanish in the limit as nf approaches 2 in the 
paramagnetic case. In fact, this effect should occur in the Kondo limit for any /" — /"^^ 
problem with n 7^ (the n = case is special since the /° state has zero f electrons). 
We suggest this as a possible reason for the lack of any definite evidence for crystal field 
splittings in UPt^. Evidence for crystal field effects have been seen in other heavy fermion 
superconductors {URu2Si2, UPd2Al^), but these are a reflection of the crystal fleld splitting 
of the p manifold, not of the quasiparticles states. We note that in these cases, the ground 
state is magnetic, causing a profound effect on the quasiparticle electronic structure, so the 
above considerations based on a paramagnetic solution would not apply (a weak moment 
state is found in UPt^, but it does not appear to affect the quasiparticle electronic structure). 

In conclusion, by studying the f^ — P Anderson lattice model, we flnd that multi-particle 
correlations have a profound impact on the quasiparticle electronic structure. The results 
of this paper should have particular bearing on uranium heavy fermion metals where such 
multi-particle ground states are known to exist. 
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FIGURES 

FIG. 1. Renormalized crystal field splitting, AA, versus nj, with a bare splitting of 0.1 eV. The 
straight line is an analytic approximation valid near an Uf of 2 derived in the text. The change in 
behavior near an n/ of 1.2 is due to a jump in the chemical potential from the lower to the upper 
5/2 quasiparticle band. 

FIG. 2. Absolute value of the separation of the renormalized f level position in the 5/2 channel 
from the chemical potential versus nj, which should be an estimate of the Kondo temperature. 

FIG. 3. f electron counts in the 3/2 and 5/2 channels versus nj. 
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